Spectral properties in the charge density wave phase of the half-filled Falicov-Kimball 
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We study the spectral properties of charge density wave (CDW) phase of the half-filled spinless 
Falicov-Kimball model within the framework of the Dynamical Mean Field Theory. We present 
detailed results for the spectral function in the CDW phase as function of temperature and U . We 
show how the proximity of the non-fermi liquid phase affects the CDW phase, and show that there 
is a region in the phase diagram where we get a CDW phase without a gap in the spectral function. 
This is a radical deviation from the mean-field prediction where the gap is proportional to the order 
parameter. 



I. INTRODUCTION 

Metzner and VoUhardti pioneered a new approach to 
the study of strongly correlated electron systems which 
is exact in the limit of infinite dimensionality. Generally 
referred to as Dynamical Mean-Field Theory (DMFT). 
the method has been developed further in the subsequent 
years and has led to substantial progress in our under- 
standing of these systems^i^. Soon after the work of Met- 
zner and Vollhard^, in a series of papers^, Brandt and 
Mielsch showed that the large dimensional limit of the 
spin-less Falicov-Kimball model (SFKM)^ is exactly sol- 
uble and studied various aspects of the solution. One of 
the first issues examined in the SFKM using the DMFT 
was that of ordering into a two sublattice, "checker- 
board", charge density wave (CDW) state that is ex- 
pected to arise in the model at half-filling on a hyper 
cubic lattice. Because the hyper cubic lattice is bipar- 
tite, a transition at a non-zero temperature is expected^ 
for all U and indeed this is true within the DMFT as 
welli. 

Since the original work^, a number of other studies of 
the 'uniform phase' at high temperatures, the CDW or- 
der and of phase separation away from half filling have 
f oUowed^i^ i 1 1 1 1 1 1 ^ ^ on both the hyper cubic and Bethc 
lattices. The homogeneous or uniform phase at high tem- 
peratures is a non-Fermi liquid^. A detailed study of 
transport in this phase has been carried out by MoUer et 
a&. Van dongcn- pointed out that for a lattice with a 
bounded density of states (DOS) the half-filled system in 
the uniform phase displays a metal-insulator transition at 
a non-zero U. Since the DMFT is in the thermodynamic 
limit, one can determine the transition temperature for 
the continuous transition from the uniform to the CDW 
phase simply by finding the temperature at which the 
charge susceptibility at wave- vector Qo = (tTjTt, ...) di- 
verges. Thus one can obtain the phase diagram of the 
half filled system purely from the knowledge of the uni- 
form solution. The phase diagram for commensurate and 
incommensurate filling (away from half filling) in one and 
two dimensions within the frame work of the DMFT has 
been studied by FreericksiS. 



Despite all this work, the properties of the CDW phase 
itself do not seem to have been explored a great deal. In 
this contribution, we study the spectral and transport 
properties of the charge density wave (CDW) phase of 
the half-filled SFKM in the DMFT approximation at a 
level of detail not reported before, to the best of our 
knowledge. We present detailed results for the spectral 
functions in the CDW phase as functions of temperature 
and U. We show how the non-Fermi liquid behavior in 
the uniform phase affects the CDW phase, and show that 
there is a region in the phase diagram, not emphasized in 
the literature before, where we get a CDW phase without 
a gap in the spectral function. This is a radical devia- 
tion from the static mean-field prediction where the gap 
is proportional to the order parameter, and the CDW 
phase is always gapped. The 'transition temperature' 
Ti{U) between this novel gapless 'phase' and the conven- 
tional gapped CDW phase has an interesting 'reentrant' 
structure as U increases. 



II. DMFT OF SFKM; FORMALISM 

The Hamiltonian for the SFKM on a lattice with sites 
labeled by i can be written as 



H 



i,j i 



It includes spin-less conduction or band electrons b with 
hopping parameters Uj , which, in this paper, we assume 
to be non zero only for nearest neighbors, and localized ^ 
electrons conserved at each site with site energy e^. /i is 
the chemical potential for the h and £ electrons. There is 
an on-site coulomb repulsion U between 6 and L In this 
paper, we confine ourselves to the particle hole symmetric 
case, with A^~^X]i^^i = ^^^Tli < '^^i and 
= N'^ili < "bj >= 1/2 where N is the 
total number of sites in the lattice. On bipartite lattices 
this is achieved by the choice = 0, /i 
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The exact solubility of the SFKM within the DMFT 
approximation arises as follows^. In infinite dimensions, 
or in finite dimensions within the DMFT. the self en- 
ergy, irreducible vertex functions etc., are purely local, 
and in diagrammatic perturbation theory, for example, 
are given by sums of skeleton graphs involving only the 
local Green's function. The problem is hence mapped 
to a single site or impurity problem embedded in a self 
consistent medium representing all the other sites of the 
lattice. In a functional integral formalism, the resulting 
single site effective action at site i is given by^ii^ 



Spff = 



Jo Jo 

+ Unu / dTh\{T)h{T) (2.2) 

Here 6J(t) and hi{T) are fluctuating fermionic Grass- 
mann fields corresponding to coherent states of the h elec- 
trons, and n^i = 1 or corresponding to the t state at 
site i being occupied or not. 5''(t) is the bare on site 
local propagator (also referred to as the host or medium 
propagator) for the self consistent single site problem. 
It is essentially the local propagator at site i excluding 
the self energy processes at that site (since they will be 
recovered by solving the single site problem) but includ- 
ing the self energy processes at all the other sites of the 
lattice in some average way. (Site indices are being kept 
track of in the equations above with a view to using them 
in contexts involving broken symmetry solutions such as 
the CDW phase.) 

Since the effective action eq. (|2.2p is quadratic in the 
variable b for a fixed na , the local Green's function can be 
calculated exactly for an arbitrary host Green's function 
as 



Gil 



Wo, 



u 



(2.3) 



Here G" = G"{iuJn) (and similarly for Q, etc.), icOn = 
i{2n + l)7rT arc the Fermionic Matsubara frequencies, 
and wii and wpi = (1 -Wi^) are the annealed probabil- 
ities for the £ state at site i being occupied and empty 
respectively, wi,;, also equal to the average i electron 
number at site i, can be calculated as 



Wli 



Oi 



(2.4) 



in terms of Zu and Zoi, the constrained partition func- 
tions obtained by summing cxp(— Se//) over all the fluc- 
tuating b configurations for fixed nu = 1 or respectively. 
It is straightforward to verify that 



exp 



(2.5) 



A renormalized or effective £ electron energy at site i 
can be defined by expressing wi; as np(eg, — fi), whence 



4, = l-i + Tln{^). 



(2.6) 



The local Green's function host Green's function 
t/" and self energy EJj' are related by the local Dyson 
equation. 



(2.7) 



Using this equation in reverse, in the form, 

(^:r'^j^ + {G:i)-\ (2.8) 



substituting into eq. (|2.3p . and cross-multiplying, one 
gets^ a quadratic equation which can be solved to ex- 
press the self energy directly in terms of the local Green's 
function and wu. The result is 

= I - ^ (1 - v/l + (C/G^O' + 4(AwiOC/Gjj) , 

(2.9) 

where we have introduced 

Aw, = (wi, - i) = (1 - wo,). (2.10) 

(The sign of the root chosen in writing eq. (|2.9p is so as 
to reproduce the known analyticity properties of the self 
energy.) In terms of the ratio Xi = Zu/Zoi which can be 
calculated, using eq. (|2.5p . as 



Xi = exp 



+ ^ In 



we can write 



Aw,- = 



1 / X,,- 1 



2 \xi + l 



(2.11) 



(2.12) 



We note that, although according to eq. (|2.9[) S" seems 
to depend explicitly only on at the same Matsub- 
ara frequency, implicitly, via its dependence on Awi and 
hence on Xi, it is in general a functional of G^ at all 
frequencies iuim- 

To complete the DMFT procedure, one has to next 
invoke the DMFT self consistency relation which deter- 
mines the local Green's functions in terms of the local self 
energies via lattice propagators, based on treating the lo- 
cal self energy as a good approximation (which is exact 
in infinite dimensions) to the lattice self energy. For this 
purpose we need to use the lattice Dyson equation. 



^Lu^ + ^,-J:^i)Gli-J2t^lG'il^6,, 



(2.13) 



The properties of this relation depend very much on the 
phase one is interested in. 
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III. THE UNIFORM PHASE AND THE CDW 
INSTABILITY 

At high temperatures the b electrons and i electrons 
are uniformly distributed on the lattice for all electron 
concentrations; there is no long range order, and no bro- 
ken symmetry. In this case of a homogeneous or uniform 
phase, wij=wi, E" = S„, etc.; i.e., the same for all the 
sites i. And at half filling, by symmetry, wi = ^, so that 
A Wi = 0. Hence I]„ becomes a single explicit function 
of Gn given by 



1 



^" 2 2G 



(l - v/1 + (UGny) . (3.1) 



The DMFT self consistency condition eq. (|2.13p is easily 
solved in this limit simply by fourier transforming. The 
lattice Green's function is 

Gk(«t^„) = (jcj„ + - I]„ - ek)~^ . (3.2) 
Hence the local Green's function is 



^o(e) 



de 



(3.3) 



where Do{e) is the bare density of states(DOS) of the 
band structure determined by iy. Most of the results 
reported in this chapter have been obtained using the 
semi-circular DOS (SDOS) characteristic of the Bethe 
lattice in infinite dimensions, for which 

Doie) - ^^4^2- £2 , |,| < 2t , (3.4) 

whence Gq{z), the complex Hilbcrt transform of the bare 
DOS, can be obtained analytically, as 



Go{Zn) 



Zn + ^Zl - 4^2 



(3.5) 



with Zn — iuJn+lJ-^'Sn- Wc also report some results using 
the actual density of states for a square lattice whence 
Go{z) has to be evaluated numerically. We will refer to 
these as 2DD0S results. In cither case, the local Green's 
function and the self energy can be obtained by self con- 
sistently solving eqns. p.ip and p.3p . 

In the uniform case at half filling the temperature com- 
pletely drops out of the problem, as the £ occupancy 
probabilities are fixed by symmetry. The metal-insulator 
transition we alluded to above^ hence takes place at a 
fixed value of U independent of temperature within the 
DMFT. For the SDOS case this value is U = 2t. The 
temperature dependence of the b electron DOS can be 
restored by including spatial fluctuations using approxi- 
mations that go beyond DMFTi^. 

As the temperature is lowered, the homogeneous (dis- 
ordered) phase becomes unstable with respect to an 



ordered phase where both the b and i charge densi- 
ties with spatial modulations at some ordering wave- 
vector q develop; i.e., ^ ■ ea;p(iq • Ri) < 6^6^ > and 

J2i exp{i(\ ■ Ri) < > acquire nonzero values. If the 
transition to this ordered phase is continuous (second- 
order), then it occurs when the charge-susceptibility (at 
the ordering wave vector) diverges. The static, or zero 
(external) frequency 6-chargc-susceptibility at the order- 
ing wave vector q is obtainable as the Fourier trans- 
form of the 6-charge-density - 5-charge-density correla- 
tion function: 



X(q,r) ^ -l^e'<i-(H,-H. 



dT <b]{T)bj{T)bl{Q)bk{Q)> , (3.6) 

oc 



n— — oo 



Here n indexes the Matsubara frequency of the fcrmionic 
propagators that arise in a diagrammatic perturbation 
theory evaluation of x- Using the generalized Dyson's 
equation^ valid for two particle propagators, one can re- 
late Xn to the simplest "bubble graph susceptibility" 
in the form 

Xn(q)=X^(q)-X°(q) £ ^^^xM ^ (3-7) 



where x^ is given by^ 

X°(q) = -^G„(k + q)G„(k), 



de- 



Gq 



(3.8) 



(3.9) 



Here Do{e) is the bare DOS as before, and all of the 
wave-vector dependence occurs via the quantity 



1 



and we use 



Go(z)= / deDoie)/{z-e) 



(3.10) 



(3.11) 



to denote the complex Hilbcrt Transform of the bare DOS 
as before. Similarly, the £-£ charge correlation function 
defined by 



(q, T)^^J2 e^'*-^^^-^") < n,,n,, > (3.12) 
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IV. SPECTRAL AND TRANSPORT 
PROPERTIES OF THE 
CHARGE-DENSITY- WAVE PHASE 

The checkerboard charge-density-wavc(CDW) phase 
has two incquivalcnt sublatticcs, A and B, defined by: 

exp{iQ ■ = +1 for i belonging to A (4.1) 
exp{iQ ■ Ri) = —1 for i befonging to B 

with Q = (tTjTTjTt, ). The two sublattices sponta- 
neously develop difi'erent electronic densities relative to 
each other, but all sites of a given sublattice are equiv- 
alent. All the quantities such as wii, S", G", etc., in- 
troduced earlier in the DMFT formalism section now be- 
come double valued with respect to i. Hence we can write 



FIG. 1: (Color online) Conventional phase diagram of the half- 
filled SFKM in the r - [/ plane. Phase A, C, and D label 
the CDW phase, the uniform non-Fermi liquid metal, and the 
uniform insulator respectively. 



wid if i belongs to sublattice A 
wib if i belongs to sublattice B 



(4.2) 



Similarly, for the local self-energy and Green's function, 
one has 



if i belongs to sublattice A 
Yj^ if i belongs to sublattice B 



(4.3) 



can be shown^ to diverge at the same temperature [for 
the same value of -^(q)] as the fo-charge-susceptibility 
discussed above. Therefore, a divergence of the 6-charge- 
susccptibility [x(q, T) in eq. (|3.7|) ] also signals an insta- 
bility towards the generation of long-range order in the £ 
charge density corresponding to the same parameter X. 



The mapping q Ar(q) is a many-to-one mapping 
that determines an equivalence class of wave vectors in 
the Brillouin zone. In the d ^ oo limit, "generic" wave 
vectors are all mapped io X = 0, since cosqj can be 
thought of as a random number between -1 and 1 for 
generic points in the Brillouin zone. All the possible val- 
ues of AT (— 1 < X < 1) can be sampled, however, using 
a wave vector that lies on the diagonal of the first Bril- 
louin zone extending from the zone center (A = 1) to 
the zone corner (A = —1). An ordering at the zone 
corner [A = — 1,Q = (tt, tt, ...)] corresponds to a two 
sublattice (checkerboard) state with the i electrons pref- 
erentially occupying one sublattice [exp{iQ,.Rj) = 1] and 
the b electrons preferentially occupying the other sub- 
lattice [exp{iQ,.Rj) = —1]. Here we only consider the 
A = — 1 case. The susceptibility x(q, ?") then diverges 
at the temperature Tc (A = — 1) and Tc as a function of U 
for the SDOS case is shown in Fig. [1] Phase C denotes 
the non-Fermi liquid metal (9S(a; = 0) 7^ and tem- 
perature independent), D is the uniform insulator and 
A is the two-sublattice or checkerboard charge-density- 
wave(CDW) phase^^. 



^ii _ / if i belongs to sublattice A /aa^ 
" " \ G^f if i belongs to sublattice B ' ^ ' 

At half filling, by symmetry, 

(wiQ - ^) = -(wi6 ^ ^) = (4-5) 
Hence, the relation eq. (|2.9p becomes 

^ = (1 - Vl + (C/G,f^)2±4(Aw)C/G,f^) , 

(4.6) 

To complete the DMFT self consistency loop, we have 
to solve the lattice Dyson equation eq. (|2.13p in the 
present, two sublattice CDW context. For this purpose it 
is convenient to introduce a two-component spinor field 

with 

bkA = X! ^» exp(ik.Rj) , 

icA 

etc., where k now only takes values within the reduced 
Brillouin zone corresponding to the A sublattice alone. 
In terms of the matrix Green's function defined in the 
usual way as 

Gk„ =« V'k;V^k » ' (4-8) 
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the lattice Dyson equation eq. (|2.13|) reduces in the 
present context to the simple 2x2 matrix equation 



where 



[{iujn + A^)! - Mk„]Gk„ = 1, 



Hence Gkn is given by 



(4.9) 



(4.10) 



(4.11) 



with 



A,B 



iuJn + /i — E 



A.B 



The local Green's functions G^'^ are just the diagonal 
components of the matrix: 



G„ — y ^ Gkn • 
k 

Hence we get, finally, 

'■B.A f Do{e)de _Go{zn) ^.A 



(4.12) 



(-iA,B ^ pB.A f Dpy^ 
" ?n / fiApB 



Sn 



(4.13) 



where Do{e) is bare DOS of the b band as before, 



Sn Sn 7 



and Go{z) is the complex Hilbert Transform of the bare 
DOS as defined before in eq. p.lip . Again, for the case 
of the semi-circular DOS the calculations become simpler 
as Gq{z) is analytically known (cf eq. p.Sp ). 

We note that Aw can be thought of as the order 
parameter characterizing the CDW phase. In the uni- 
form phase, Aw= 0, and = S^, etc., and all the 
above equations reduce to the corresponding equations 
discussed earlier for the uniform case. 



V. RESULTS AND DISCUSSION 

Next we discuss in detail our results for the DMFT 
spectral and transport properties in the CDW phase of 
the half-filled SFKM. For convenience we measure all en- 
ergies such as Lo, U , etc. in units of t for the SDOS, and 
2t for the 2DD0S. Consequently, the bare h band width 
on the Bethe lattice or on the square lattice is = 4. 

The main result from our study is summarized in the 
new 'phase diagram' for the half-filled SFKM, shown in 
Fig. O As noted before in the introductory section, there 
is a region in the phase diagram, marked B in Fig. [2] and 
not emphasized in the literature before, where the CDW 
phase order parameter is non-zero, but there is no gap 
in the spectrum. This is unlike the Hartree-Fock (HF) 
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FIG. 2: (Color online) New phase diagram of the half-filled 
SFKM in the T - U plane. Phase A, B, C, and D represent 
gapped CDW, gapless CDW, uniform non-Fermi liquid metal, 
and uniform insulator respectively. 
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FIG. 3: (Color online) Variation of the order parameter (Aw) 
with temperature for difi'erent values of U. (All the energies 
are scaled in units of t.) 



mean-field theory prediction, where the CDW phase is 
always gapped. We will refer to this as the "gapless CDW 
phase". Furthermore, the 'transition temperature' Ti(U) 
between this novel 'phase' and the conventional gapped 
CDW phase has an interesting reentrant structure as U 
increases, as depicted in Fig. [2l 

This emerges from a study of results such as the ones 
shown in Fig.s [3]l5l Fig. [3] shows the CDW order pa- 
rameter (Aw ) as a function of the temperature T for 
five values of U = 0.5, 1.0, 1.5, 2.0, 2.5. In Figs. H and E 
we have shown the evolution of the local spectral func- 
tion for SDOS and 2DD0Si^ respectively as a function 
of temperature and U. In the uniform phase, the spec- 
tral functions are temperature independent as we have 




FIG. 4: (Color online) Diagrams showing the evolution of local spectral function for SDOS as a function of temperature for 
three values oi U — 0.5, 1.0, 1.5 in the first, second, and third columns respectively, and in increasing order of temperature. 
Dotted lines are for sublattice A and long dashed lines are for sublattice B. 
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FIG. 5: (Color online) Diagrams showing the evolution of local spectral function for 2DD0S as function of temperature for 
three values of f/ = 0.5, 1.0, 1.5 in the first, second, and third columns respectively, and in increasing order of temperature. 
Dotted lines are for sublattice A and long dashed lines are for sublattice B. 
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remarked before, and their U dependence can be seen in 
the bottom panels of Fig. [4] for the SDOS case and in 
Fig. for 2DD0Si^. Basically, as U increases, there is 
a transfer of spectral weight from low frequencies out to 
large frequencies, leading to the formation of 'Hubbard 
bands' and a 'pseudo gap'. Eventually, for U > Umi 
{Umi = 2 for SDOS) a real gap opens up, and a metal 
insulator transition ensues, as indicated in Fig. [21 As we 
decrease the temperature at a fixed U within the checker- 
board phase, there is spectral weight transfer from low 
frequencies to near a threshold frequency corresponding 
to the HF gap. However, as can be seen in the next (sec- 
ond from bottom) set of panels in Figs. [4] and [5l a real 
gap in the spectral function does not develop just after 
one crosses into the checker-board phase as would be pre- 
dicted, for example, by the HF approximation. Rather it 
develops only below the lower critical temperature Ti{U) 
shown in Fig. [5) Even more interestingly, Ti{U) has a 
reentrant feature as function of U (see Fig. . 

We note that at T = 0, Aw=0.5, = U and 0. 
The DMFT spectral functions at T = are thus identi- 
cal to the HF spectral functions, corresponding to split 
bands with dispersion ±-\/ {U /2)2 -|- e^. The gap in the 
spectral function is hence exactly equal to U. However, 
as we increase the temperature the value of Aw devi- 
ates slightly from 0.5 and the spectral function for each 
sublattices develops two peaks around a; = and one 
observes three gaps as can be seen in the topmost pan- 
els in Figs, m and [5l As we increase the temperature 
further, these two peaks grow in intensity and the gaps 
away from to = close up and only the gap around uj ~ 
remains. Eventually, past the lower critical temperature 
Ti{U) mentioned above, even the gap around uj = van- 
ishes, but the order parameter is still non-zero, until T 
crosses Tc whence Aw becomes zero and the local spec- 
tral function becomes uniform. This behavior is generic 
for U < Umi- For U > Umi, the spectrum is always 
gapped. 

In the gapless CDW phase B, the b electrons are in a 
state which allows gapless excitations and coexist with £ 
electrons in a gapped state. We can view the effective 
£ electron energy as the centroid of the i electron spec- 
tral function. The gap in £ electrons spectral function is 
therefore definable as e^^ — . The absence of a gap 
in the 6-electron spectral function arises from the same 
source that causes the 6-electrons in the uniform phase 
to acquire a non-Fermi liquid character, namely the £- 
electrons which act as a source of disorder scattering, 
leading to a finite life time of the b electrons at the Fermi 
energy. In the CDW phase close to Tc and for U not 
too large, although there is long range staggered order in 
rii, there are large thermal fluctuations in ni which again 
act as disorder scatterers for the 6-electrons, leading to 
a finite life-time at the Fermi surface and to a gapless 
6-electron spectral function. This feature goes away at 
low temperatures or for large U, as the order in ni gets 
stronger and the fluctuations get reduced. 

To get a quantitative estimate of the lower critical tem- 
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FIG. 6: (Color online) Variation of r(= 3 (T.a + 'Sb) at lo =0) 
with temperature for different values of U for SDOS. 
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FIG. 7: (Color online) Variation of F with the order parameter 
(Aw) for different values of U for SDOS. 

perature T;, we have studied the quantity 

F = 3 (Sa + ^b)lj=o 

as a function of temperature, as the vanishing of F sig- 
nals the opening up of the gap around uj = 0. Figs. [H 
and [7] depict F as a function of temperature and Aw re- 
spectively for U = 0.5,1.5,1.55,1.75. We note that at 
U = 0.5 the gap closes up at T/ = 0.029 (corresponding 
Aw=0.240). As we increase U the gap closes up at higher 
values of temperature, i.e., at T/ = 0.052 (Aw=0.435) 
for [7 = 1. However, as we increase U further, the gap 
closes up at lower values of Ti up to U very close to a/2 
and Ti = 0.0108 (Aw~ 0.5). Beyond [/ ~ \/2 the lower 
critical temperature T; for the closing of the gap again in- 
creases to join the versus U curves at [/ = 2.0 . Hence 
we get the reentrant curve for Tj as shown in Fig. [21 
The lower critical temperature Ti shown in Fig. [2| was 
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FIG. 8: Variation of the critical temperature Ti with U. 



FIG. 9: Variation of the order parameter (Aw) with the crit- 
ical temperature T;. 



obtained by analyzing the data of spectral functions and 
the order parameter as discussed above. Using the rela- 
tion Y.^{uj = 0) = —{Ti^{uj = 0))* and the expression for 
Yi^iuj = 0) from eq. (j4.6p . and invoking the condition 
for the vanishing of F, we obtain the following implicit 
equation for Tf. 



^ATi) + ^Re{G^{oj = 0; T,)) = 0. (5.1) 

Ti(U) can be very accurately determined using eq. (j5.ip 
and the results are exhibited in FiglH and in the inset of 
Figl2] We have also computed the order parameter (Aw) 
as a function of T; and U on the curve Ti{U), and we have 
shown these in Figsin]and[TUl As we move along the curve 
Ti{U) in Fig. [2] the order parameter (Aw) first increases, 
helped both by the increase of U and the decrease of 
T, approaching a maximum value very close to 0.5 at 
U~ a/2- However, as we move to larger values of U, the 
order parameter (Aw) starts decreasing and vanishes at 
U=2.0, where T/ joins T^. 

When Aw is very close to 0.5, then and are 
very close to U and respectively. Substituting these 
values in eqs. (|4.13p and (|5.1|1 . we obtain 



u 

2 



4t2 



-U. 



(5.2) 



Solving eq. (|5.2p . we can verify that the dip of the 
lower critical temperature Ti occurs 

It is interesting to ask how these novel features in 
the spectral functions affect experimentally measurable 
quantities. In particular, we have studied the optical 
conductivity a-{u!) as a function of T and U. In the CDW 
phase it can be shown to be given byi^ 



.(.) |^d.A,.(.)y_^T.(a^A(.,.' 



Aw 



0.5 



0.4 



0.3 



0.2 



0.1 




0.5 
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u 



1.5 



FIG. 10: Variation of the order parameter (Aw) with U on 
the curve Tj. 



„. , , ..np(uj') — nn.(iLj' + uj) 
xa'^Aie^Lo' + uj))-!^ — -, 



(5.3) 



where is the Pauli matrix, A(e,w') 



1 c> 



3G,(t^'), 



is the matrix spectral function for the matrix Green's 
fimction given by eq. (|4.1ip . and Dtr{e) is the transport 
DOSi^iii. We plot (y{u!) for different values of temper- 
atures for a fixed value of U = 0.5 in Figs. [TT] and [12] 
for the SDOS and 2DD0S respectively. When T=.016 
one can see from the spectral function (first diagram of 
column 1 of Fig. |4]) that the gap around w = is roughly 
0.25, and there are two bands of low spectral weights 
which are separated from the main bands by a second 
smaller gap (0.05) around uj = 0.5. These features are 
reflected in a{iu) (see Fig. [11] dot-dashed curve) which 
shows that there is no optical response up to a; = 0.25. 
Then there is a rise, a small dip around lu = 0.5, fol- 
lowed by a sharp peak. As we increase the temperature 
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FIG. 11: (Color online) Optical Conductivity cr{ij) at different 
temperatures for a fixed value of (7 = 0.5 for SDOS 
100 




--- T=0.016 

— T=0.026 

— T=0.030 

— T=0.034 



^ 



0.125 0.25 0.375 0.5 0.625 0.75 0.875 1 

(0 



FIG. 12: (Color online) Optical Conductivity cr{uj) at different 
temperatures for a fixed value of i7 = 0.5 for 2DDOS 



to T = 0.026 (see second diagram of column 2 of Fig. d]), 
the spectral function now has a very small gap around 
u! ~ 0, this feature of spectral function is reflected in the 
corresponding a{u!) (see Figs. [TlT long dashed curve) and 
[T^ (dot-dashed curve)). Similarly, the other two curves 
for the optical conductivity at T = 0.03 and T = 0.032 
(Fig.[TT]dottcd and solid lines respectively) correspond to 
the spectral functions in the last two diagrams of column 
1 in Fig. 31 We note however that the dc conductivity 
a{uj = 0) does not capture these features, and reflects 
only the transition from the CDW to the uniform phase 
which shows up as a slope discontinuity at T^., as can be 
seen in Fig. [T51 



FIG. 13: (Color online) dc-conductivity as a function of tem- 
perature for a fixed value oi U = 0.5, 1.0, 1.5. 



VI. CONCLUDING DISCUSSION 

In conclusion, we have presented results from a de- 
tailed DMFT study of the spectral functions in the CDW 
phase of the half-filled SFKM as function of tempera- 
ture and U. We have shown that the proximity of the 
non-Fermi liquid metallic phase afi^ects the CDW phase, 
leading to a region in the phase diagram where we a get 
CDW phase without a gap in the spectral function. In- 
terestingly, this gapless CDW phase shows a reentrant 
transition to the gapped CDW phase as U increases. 
This is a radical deviation from mean-field prediction 
where the CDW phase is always gapped, with the gap 
being proportional to the order parameter. We have also 
discussed how these features affect response functions, 
e.g., the optical conductivity. It would be interesting to 
study whether, and to what extent these features sur- 
vive when one goes beyond DMFT for the SFKM, e.g. 
in more sophisticated approximations such as the Dy- 
namical Cluster Approximationi^ii^, clustcr-DMFT— or 
Variational Cluster Approximation^^, which include the 
effects of short range inter-site correlations. 

We note that gapless CDW phases are easy to achieve 
even within the Hartree approximation by considering 
second neighbor hopping, which gets rid of nesting. What 
is novel about the gapless CDW phase discussed here 
is that it is correlation induced, and appears despite 
the presence of perfect nesting. It is interesting to ask 
whether such phases can appear in other models with 
strong correlations. The normal (repulsive) Hubbard 
model does not have CDW instabilities. But an extended 
Hubbard model with nearest neighbor repulsion V would. 
Wc believe that at intermediate values of the Hubbard U 
and appropriate values of V the extended Hubbard model 
with nearest neighbor hopping could exhibit a gapless 
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commensurate CDW phase. For, the DMFT treatment 
of such a model in the presence of a CDW would cor- 
respond closely with the recent study^i of correlation ef- 
fects in a two sublattice band insulator where correlations 
were shown to induce metallicity. 

On the experimental front, gapless CDW phases have 
recently been observed in 2H- transition metal dichalco- 
genides (eg., ref. [l^ and references therein). In these 
compounds, strongly coupled electronic and lattice de- 
grees of freedom are involved in the generation of the 
CDW instability, and the gapless feature has been at- 
tributed to very large second neighbor hopping. Since in 
our model we are considering only electronic degrees of 
freedom, and only nearest neighbor hopping, we have not 
compared the experiments directly with our findings. 
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